EXISTENCE OF RADIAL STATIONARY SOLUTIONS FOR A 
SYSTEM IN COMBUSTION THEORY 



(1.1) 
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Abstract. In this paper, we construct radially symmetric solutions of a non- 
linear noncooperative elliptic system derived from a model for flame balls with 
radiation losses. This model is based on a one step kinetic reaction and our 
system is obtained by approximating the standard Arrehnius law by an ig- 
nition nonlinearity, and by simplifying the term that models radiation. We 
prove the existence of 2 solutions using degree theory. 



1. Introduction 
This paper deals with radial solutions of the system of equations 

An - en = -vf(u) in R 3 
Av = vf(u) in R 3 
u > 0, v > in R 3 



lim u(x) = 0, lim v(x) = 1, 



where e > and / is an ignition type nonlinearity, that is, there exists < 6 < 1 
such that 

(i.2) f(t) = o vt < e, f(t) > o vt > e. 

System (jl.ip arises as a model problem for some reaction diffusion systems in 
combustion theory. We will describe the connection with these models at the end 
of the section. 

We call (u,v) = (0,1) the trivial solution of (|1.1|) . Note that for any solution 
(u, v) of (jl.ip one has 

< u < 1 and < v < 1, 

as can be deduced from the maximum principle. 

Wc arc mainly interested in existence and multiplicity of nontrivial radial solu- 
tions of (jl.ip . The first observation in this direction is that for e > large (|1.1|) 
has no nontrivial solutions. Indeed, by the maximum principle u + v < 1 and then 

-Au < (1 - u)f(u) - eu in R 3 . 

If e > is large then (1 — u)f(u) — eu < for all u G [0, 1] and we conclude that 
u = 0. Better estimates for the quantity 

e* = sup{e > : (jl.ip has a nontrivial radial solution} 

are given in Section [6l 
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Suppose that u, v is a nontrivial radial solution. Then for some r > we must 
have u(r) > 9. Otherwise f(u) = and then Aw = 0, and by the Liouville theorem 
v = 1. Since < u + v < 1 this would imply u = and then u, v is trivial. So for a 
nontrivial solution u, u, since lim r _ i . 00 w(r) = 0, there is /3„ > such that 



The parameter (3 U , which we will write simply as /3, then serves to distinguish 
different solutions. 

Theorem 1.1. Assume f : [0,1] — > [0, +oo) is continuous, satisfies (|1.2I) and for 

some C > 



T7ien e* > and i/iere exist < e\ < e* such that for < e < E\ there are at least 
2 solutions of (|1.1[) . One of them has bounded /3 as e — > and the other has (3 in 
the range 5/\J~£ < (3 < l/(6' v / e) as e — > 0, where S > is fixed. 

A very natural and interesting question is the stability of the solutions con- 
structed in Theorem 11.11 Based on the works [101 E] we conjecture that the 
solution with bounded f3 is unstable and that in part of branch of solutions with 
large (3 the solution is stable, at least with respect to radial perturbations. 

In some cases one may want to consider a discontinuous nonlinearity, such as 
the Heavisde function f{u) = X[u>6]- With this example in mind we introduce the 
following hypothesis 

(1.4) / is continuous in (9, 1] and lim tl _ i .e + f(u) exists. 

Theorem 1.2. Assume f : [0,1] [0,+oo) satisfies (fTT^j) . (jOj) . and (flU) . Then 
e* > and there exist < E\ < e* such that for < s < £\ there are at least 2 
solutions of p. II) . One of them has bounded ft as e — > and the other has /3 in the 
range S/^/s < /3 < \/(9^/e) as e — > 0, where S > is fixed. 

The solution constructed in Theorem ll.2l is such that the set {r £ [0, oo) : u(r) = 
9} is finite, and hence the equation (|1.1[) holds a.e. in M 3 . 

The motivation to consider a discontinuous nonlinearity is only mathematical. 
However the example f(u) = X[u>e] is interesting since it provides a situation where 
explicit calculations are possible. Theorem 11.21 shows that in part the conclusions 
obtained for f(u) = X[u>e] remain valid for more general non-linearities. 

For f(u) = X[u>8] explicit calculations lead to an equation for f3 and e in order 
for a radial solution to exist. In Figure [T] we show the numerical solution for this 
relation when 9 = 0.5, with ft in the vertical axis and i/e in the horizontal axis. 
It shows that for < e < e* there are 2 solutions. Solutions in the lower branch 
satisfy u > in [0, (3) that is, the reaction takes place in the ball of radius (3. The 
same happens for points in the upper branch which are to the right of the special 
point marked in the graph. To the left of that point the solution satisfies u > 6 in 
an annulus of the form r € (/? — Lp, (3). Thanks to the explicit form of the relation 
bewteen /? and e we can compute the asymptotic behavior of the curve as (3 — > oo, 
and we find that 



U (I3 U ) = 9 and u{r) <9 Vr > A 



(1.3) 



f(t)<cf(u) wt<u, t,ue[e,i 




lim \fif3 — ao > and lim Lp = Lq > 0, 
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Figure 1. Bifurcation diagram for / = X[u>6] with 9 = 0.5; j3 is 
in the vertical axis and \fe in the horizontal axis 



where ao,£o is the unique solution of the system of equations 



Because of the information on the Heaviside nonlinearity one can conjecture that 
for general / there should be a similar relation for j3 and e as j3 — > +oo. We present 
in Section [2] nonexistence results for general ignition nonlincarities satisfying (|1.2j) 
and (jl.3p , that capture this relation, and roughly speaking say that no solution can 
exist if y/efi is either too large or too small, provided /3 is taken large enough. Using 
these nonexistence results and degree theory we can give a proof of Theorcm ll.il 
This is done in Section [3J In Section @] we give the proof of Theorem 11.21 by 
approximating the discontinuous nonlinearity by continuous ones. Section [5] is 
devoted to the explicit computations for the Heaviside function. Finally Section [5] 
contains a finer estimate of e*. 

As mentioned before, system (jl.ll) arises in connection with some models in 
combustion theory, more precisely, in the flame ball problem for a weakly premixed 
gas sensitive to radiative heat losses. In such a mixture, it is known that apparently 
stationary spherical structures appear, which are called flame balls [271 1501 [TU1 [TI] . 
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In [27J [301 HOJ HH EU E21 [29] the following reaction diffusion system has been 
used to model a combustion process where flame balls can appear: 



(1.5) 



P C P — = V • (AVT) 



dT 

~di 

dY 



T 



QBYp __b_ 
-e RT 



- BYpe 
T 

-L OO 



q(T) 



as 



where T is the temperature, Y the reactant concentration, Y^ > 0, > arc 
the reactant concentration and tcmpereature at infinity, and C p , R, Q and m are 
respectively the specific heat capacity at constant pressure, the perfect gas constant, 
the chemical heat release and the molecular mass of the reactant. The term q(T) 
represents radiative losses. The reaction is characterized by the one-step Arrhenius 
kinetics ~ Be^r where B is a constant. Furthermore, the hydrodynamics effects are 
neglected, i.e the density p, the thermal conductivity A and the diffusion coefficient 
p are constant. See also numerical simulations in [7] I21 [ [22 ] . 

After the seminal work [34], the traveling front problem for systems like (|1.5j) 
has been investigated by several authors, for example, [16 ] [17 ] 131 ] . 

In absence of radiation, i.e. when q = 0, there are many works dealing with 
(|1.5[) . see for instance [3] [6] [25] [28] [32] and the references therein. Also we remark 
that the stationary version of (|1.5[) without radiation, leads to system (|1.1[) with 
e = which reduces to a scalar equation, since v = 1 — u. There is a huge amount 
of literature concerning existence of radial ground states for semilinear equations, 
so we mention here only some classical references [5] I24| on the problem in entire 
space. When the problem is treated in a bounded domain see the book [T] and 
|15[ I23j for multiplicity results in the case of Arrehnius non-linearity. The paper [9] 
contains interesting numerical computations of the bifurcation diagram in the case 
of the full coupled system in a interior of a sphere. 

A common simplification of (|1.5[) under the assumption of large activation energy, 
that is E » 1, is to assume that the source term for the reaction is concentrated 
on a very thin layer, typically a sphere. This approach is taken for example in 
[TU1 im 121)1 12"9"] and leads to the free boundary problem 



(1.6) 



' dT 
~di ~ 


AT 


dY 


1 


~di ~ 


L 

Le 


Y = 





q 







Y 

1 OO 



' -Be~ — 5{r-R(t)) 
B(0,R(t)) 
E 3 \ 5(0, R(t)) 
T^T^ as Id -)■ 



where R{t) is the radius of the front where the reaction takes place, 5 is the Dirac 
measure and T* is the front temperature. In (TU] [TTJ [T2] the authors analyze 
the stability of stationary solutions of (|1.6[) . In a similar framework, existence 
and stability of flame balls and travelling flame balls have also been studied in 

Ha EH QU EES 0- 

We arrive at (|1.1[) by introducing the following simplifications: 
a) assume the radiative loss to be linear, i.e., q(T) = a(T — T^) where a > 0, 
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b) approximate e «r by r)(T — To)BYe St where To is an activation tem- 
perature and 77 is a cut-off function satisfying i] > in R + and rj = in 



As in [551 [33] one can model radiative heat losses using Stefan's law q(T) = e(T 4 — 
T^) for some constant e > 0. When T is close to we can write q(T) « 
AeT^iT — Taa). As a step towars understanding more general situations, we assume 
that this linear relation holds for all T, that is, we assume a). Other linear or 
piecewise linear approximations have been used before, for instance in 8 , 29 [ I33 j . 
Assumption b) corresponds the a standard approximation in combustion theory to 
avoid the cold boundary difficulty, see [4]. 

After introducing dimensionless variables u, v, corresponding to temperature 
and reactant concentration, the stationary version of ([1.5)1 becomes 



(1.7) 



' Au + vg{u) — cu = 
- vg{u) = i 
0, v —> Voo as 



1 . 
— Av 

Le 



where Le > is the Lewis number, c > 0, > and g is an ignition type function, 
that is, there is 9 > such that g(u) — if u < 8. 

We stress that our results are valid for any value of Le > 0. Indeed, since we are 
considering stationary solutions, the following change of variables will allow us to 
assume that Le = 1. Letting v(x) = v(^/Lex), u(x) = ^u(-\fLe~ x) transforms 
system (jl.7[) into (jl.lj) , where e — > and f(u) = g(j^u). We observe that / 
is still is an ignition type nonlinearity. 



2. Apriori estimates 

The purpose in this section is to establish nonexistence results in some ranges of 
the parameters. 

Given a nontrivial solution u, v of (|1.1[) let /3 > be such that 
u(fi) = 9 and u(r) <6 Vr > f3. 
Setting u(r) = u(/3r), v(r) = v(f3r) these new functions satisfy 



(2.1) 



d ( 2 du i 
or \ or 



-r 2 fi 2 vf{u) 



d 



,dv 



r 2 /3 2 vf(u) in R+ 

lim u(r) = 0, lim v(r) = 1 

- r— >oo r— >oo 

In the sequel we will study ([2.1)1 in the following set of functions 

S = {(u, v) : u, v are C 1 ([0, 00)) and u(l) = 8 and u(r) < 9 for r > 1}. 
Let ho '■ K — > R be such that 

(2.2) h o (t)=0 Vt<9 and h (t) > Vt > 9 

(2.3) ho is nondecreasing. 
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We consider now functions / : R — > R such that 

(2.4) f(t) = Vt<9, 

(2.5) f>ho in R 

(2.6) /(i) < C /(u) V0 < t < u < 1. 

Lemma 2.1. There is no solution in S to (|2.ip «/ 
Proof. Since 

Au - e(3 2 u = for r > 1, = 0, Um u(r) = 0, 

r— >+oo 

we have an explicit formula 

e -(r-i)/3V£ 

it(r) = 9 for all r > 1. 

r 

^From this we find 

(2.7) tt'(l) = -0(l + Py/e). 
Similarly, since Av = for r > 1 and lim r -s.oo w(r) = 1 we have 

T 

v(r) = 1 Vr > 1 

r 

where < 7 < 1. This yields 

i/(l)=76 (0,1). 
Integrating the equation for w in (0,1) implies 

(2.8) (3 2 [ s 2 vf(u)ds = v'(l)=~f. 

Jo 

Let ro £ [0, 1) be such that u'(ro) = 0. Then integrating the equation for u in 
(r ,l) yields 



u'{l)=e0 2 s 2 uds~l3 2 s 2 vf{u)ds 



s 2 uds~l3 2 I - 2 ~ 

To Jro 

This formula together with (|2.7p and (|2.8p gives 

v'{l)>!3 2 ( s 2 vf(u)ds = ef3 2 f s 2 uds + 9{l + 

J r Jr a 

and it follows that 

< 1. 

□ 

Lemma 2.2. There is /3q > 0, S > depending only on ho, 9, Co such that for all 
ft > /3o, — ^ arl ^ a ^ / satisfying (|2.4[) , (|2.5p . (|2.6p £/iere is no solution in S to 
the system (|2.ip . 
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Proof. We treat the case e > since the situation e = is similar. 
As before 

e -(r-l)/3v^ 

u(r) = 9 for all r > 1. 

r 

^From this we find 

(2.9) tt'(l) = -6(l+j3y/e). 
Similarly 

v(r) = 1 Vr > 1 

r 

where < 7 < 1. This yields 

i/(l)= 7 6(0,l). 
Integrating the equation for u we see that 

(2.10) r 2 u'(r) =e/3 2 [ s 2 uds-(3 2 [ s 2 vf(u)ds. 

Jo Jo 

^From this and (|2.9p it follows that 

(2.11) /3 2 [ s 2 vf{u) ds = 9(1 + psfe) + e/3 2 [ s 2 uds < 6 + 5(6 + 6). 

Jo Jo 

Integrating the equation for v in (0,1) implies 

(2.12) f3 2 [ s 2 vf(u) ds = v'(l) = 7 

Jo 

and combined with (|2.11l) yields 

(2.13) j < 9 + 5(9 + 5). 
Integrating the equation for v in (0,r) we obtain 

(2.14) v '(r)=/3 2 [ (s/r) 2 v(s)f(u(s))ds. 



Integrating this on (0, R) with < R < 1 yields 



v(R) - v(0) = (3 2 [ [ (s/r) 2 v(s)f(u(s))dsdr 
Jo Jo 

r~ 2 drs 2 v(s)f(u(s))ds 

f R 

= P 2 / (s~ l - R- 1 )s 2 v(s)f(u(s)) ds. 
Jo 

In particular, with R = 1 

(3 2 [ (s- s 2 )v(s)f(u(sj) ds = v(l) - v(Q) < 1 



This and (|2"7T2"|) give 



f3 2 / sv(s)f(u(s))ds < 2. 
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Now going back to (|2.14j) we obtain 
(2.15) 

v'(r) = f3 2 ( {s/r) 2 v(s)f{u{s)) ds < 2 I (s/r)v(s)f(u(s)) ds < - < r < 1. 



Now consider the function z = u + v which satisfies Az = e/3 2 u in M 3 . Integrating 
the relation (s 2 z')' = sf3 2 s 2 u in (0, r) 

z'(r) = e(3 2 [ (s/r) 2 u(s)ds 
Jo 

so that 

(2.16) < z'(r) <S 2 r r> 0. 
It follows that 

(2.17) -- < u'(r) < S 2 r < r < 1. 

r 

Observe that integrating (|2.16[) on (r, 1) we find 

< z(l) - z(r) < S 2 

and since z(i) = + 1 — 7 

1 + 8--/-5 2 < z(r) < l + 6» -7 0<r<l. 
In particular, using f)2 . 13[) . 

(2.18) v = z - u > 1 + 9 - 7 - S 2 - u > 1 - 8(6 + 28) - u < r < 1. 

Step 1. For any r > there exists ft(r ) depending on r , /lo an d 6 only such 
that 

(2.19) maxa<^- > ft. 

[ro,l] 2 

To prove this, suppose that max[ rn j] w > 3±i an d let r € [foil] be such that 

from (|2"T7|) 

|u'(s)|<M r < s < 1. 



it(r) = i±2. Let M = 2/r so that from (|2"T7|) 



Then 
(2.20) 

l + 6> l-6» . . 1 + 6* 1 — 6> , , l-6» 

— — < u{s) < — ^- + Vs e [r 0) 1] such that |s - r\ < 

Using (|2.12[) and since / > ho 

l>/3 2 f s 2 v(s)f(u(s))ds> [3 2 f s 2 v(s)h (u(s))ds. 

Jr Jro 

If r < ±±m i e t / = [ r _ i f ^ j r ] an d if r > l±na lct 7 = r + Using (12+81) 

and ([2~20l) : 

l>/3 2 [ s 2 v{s)f(u(s))ds 
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We take S > small only depending on 9 such that 

(2.21) i_^ + 2(5 )_^>i_f 
Then 

n9 ,(l-6>) 2 , (1 + 6 1-9 
and the claim follows. 

Step 2. For any ro > there exists C > such that 

The constant C depends only on r and 9. The conclusion from this is that 

(2.22) limsupmax(u - 6) < 0, 

and this is uniform with respect to u and /. 

As before, setting M = 2/r we have \u'(s)\ < M for r < s < 1 by (|2.17p . Let 
€ [?"o, 1] be a point such that u(r m ) = max[ ro j] u and let us write u m = u(r m ) = 
max[ ro .!] u. Then 

^(r)>^^+0 Vre[r 0; l] such that \ r -r m \<^^-. 
Since / > ft, , using (|2~T2|) . (f2~18l) . (j2~T9l) . (pT2"T|) we obtain 

l>f3 2 [ s 2 v{s)f(u(s))ds> f3 2 [ s 2 v(s)h (u(s)) ds 
Jo Jo 

>^^o(^ + ^) (l-,(* + 2*)-^ 



>P'ri— -ft 



2M V 2 

Thus there exists C depending only on rp, such that for j3 > (i\ 



This proves the claim. 
Step 3. For < a < 1 



l 

(2.23) / r 2 (u') 2 dr < 2a + S 2 +max(u - 9). 

Ja Ml 

Indeed multiplying the equation 

(2.24) -(rV)' = -eP 2 r 2 u + /3 2 r 2 vf(u) 
by (u — 9) and integrating over (a, 1) we get 

i i r 1 r 1 

r 2 (u') 2 dr = rV(u-6>) - efi 2 I r 2 u{u - 9) dr + f3 2 r 2 vf(u)(u-9)dr. 

a J a J a 

By (JUT) 

r 2 u'{u -9) = -a 2 u\a){u{a) - 6) < 2a. 
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Since 11 < 1 we have 

el 

j2 I 2 



e/3 2 \ r 2 u(u ~9)dr< 8 2 

J a 

Finally using (|2~T2|) 

/3 2 / r 2 vf(u)(u — 9)dr<ma,x(u — 9)(3 2 / r 2 v f (u) dr < max(u 

Ja [M] J a Ml 



Step 4. We finish the proof of the Lemma using a modification of Pohozaev's 
identity. Let < a < 1 be a fixed number. Multiplying (|2.24[) by (r — a)u' and 
integrating over (a, 1) 

(2.25) 

i ,i ,i 

(r it')' (r — a)u' dr = e/3 2 / r 2 (r — a)uv! dr + (3 2 / r 2 (r — a)vf(u)u' dr. 

J a J a 

A computation shows that the left hand side is given by 



1 



i i f i 



(r 2 u'y(r-a)u'dr = -r 2 (r - a)(u') 2 +-/ r 2 {u') 2 dr-a r{u') 2 dr 



(2.26) 

Thanks to (|2~23| wc find 



(1 -a)u'{l) 2 + - / r 2 {u') 2 dr-a r(u'Y dr 



a r(u) 2 dr < / r 2 (u') 2 dr < 2a + 5 2 + max(u - 0). 

Since u'(l) = —9(1 + y/e/3) we obtain from (|2.26|) . (|2.23|) and the previous estimate 
(2.27) 



(rV)'(r - a)u'dr - -(1 - a)9 2 (l + ^i/3) 2 



< 4a + 2(T + 2 max(w - 

[a,l] 



Let us compute 



i 1 



e/3 2 / r 2 (r - a)uu dr = -e/3 2 r 2 (r - a)u 2 e/3 2 / (3r 2 - 2ar)w 2 dr 



Recalling that y/e/3 < 6 this shows that 



(2.28) 
Define 

Then 



e/3 2 / r 2 (r — a)uu'dr 



< 36 2 . 



F(u)= / f(t)dt. 



/3 2 / r 2 (r - a)vf(u)u dr = p 2 r 2 (r - a)vF{u) - /3 2 / (3r 2 - 2ar)vF(u) dr 

a 

"1 



(3 2 / (r d - ar 2 )v'F(u)dr 



f3 2 (3r 2 -2ar)vF(u)dr - ft 2 (r* - ar 2 )v'F (u) dr 
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These terms can be estimated using (|2.6|l : 

F(u) < c f(u)(u -e) ye < u < 1 

and hence by ([2TT2")l 

(3 2 [ r 2 vF(u) dr < C [ r 2 vf{u){u - 9) dr < C max(u - 9) 

Ja Ja Ml 

Similarly 

a/3 2 [ rvF{u) dr < 1 [ r 2 vF{u) dr < C max(u — 9). 



To estimate the remaining terms we use ([2,15jl 



l r i 1R rl 



j3 2 r 6 v'F{u)dr <2f3 2 r 2 F{u)dr< -/ r 2 vF{u)dr 

J a J a 1 ~ V J a 

where the last inequality is consequence of (|2.18[) . f|2 . 19[) and (|2.21l) provided j3 > 
fli (a) (a £ (0, 1) is fixed). Hence 

/3 2 / r 3 v'F(u) dr < max(u - 9), 

Ja 1-9 [a,l] 



and this implies also 

2 I ar 2 v'F(u) dr < /3 2 / r 3 v'F(u) dr < 1 max(w 

1 — 9 [a,l] 



Therefore 

„2 



(2.29) 



P r (r — a)vf(u)u' dr 



< K\ max(u — 

[a,l] 



where K 1 = 5C + f^f. 

iFrom d22U), P^Tjl . (|2~2"gj) . (|2~2U)) we arrive at 

(2.30) l{l-a)9 2 (l + ^/e(3) 2 <4a + 55 2 + (K 1 +2)m&x(u-9) V/3>/3i(a). 

2 [a,l] 

We fix 5 > even smaller if necessary so that 

55 < 77 ~~r ■ 

~ 3 4 



Then fix < a < | such that 



Then for (j23fj|) yields 



10 2 

4a < . 

" 3 4 



777- < (A'i + 2) max( U -9) V/3 > ft (a) 

O 4 [a,l] 

which is not possible for (3 large enough by ()2.22j) . □ 
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3. Existence of solutions 

In this section we prove Theorem 11.11 and throughout it we assume that / is 
continuous and satisfies (|1.2p and (|1.3p . 
We work in the Banach space 

X = {(u, p);ue C lia ([0, 1]), u(l) = 0, p e R} 

endowed with its natural norm 

\\(u,P)\\ x = \\u\\c^ im) + \P\, 

where < a < 1 is fixed. 
Define for t e [0, 1] 

ft(u) = t.f(u) + (1 - t)(u -0)+ t€ [0, 1], 

where s + = max(s, 0). Let us consider (|2.1[) with nonlinearity ft that is 

{Au - eP 2 u = -p 2 vf t (u) in M 3 
Av = p 2 vf t (u) inM 3 
lim u(x) = 0, lim v(x) = 1 

To apply the non-existence results of the previous section we need to exhibit a 
function ho satisfying (|2.2p and (|2.3p . For this purpose define 

g(t) = inf {r G [0, 1] : / > t in [r, 1]} for all i > 0, 

and 

h(r) = inf {i > : ff(t) > r} for all r £ [0, 1]. 
The following properties then follow from these definitions: 

a) g is strictly increasing, continuous from the left and satisfies lim t _j. + g(t) = 

e. 

b) h(r) = for r G [0,0], h(r) > for r G (0,1], ft is nondecreasing and 
continuous in [0, 1]. 

c) / > ft in [0,1]. 

The function 

/io = min(/i, (u — 0) + ) 
satisfies and and the nonlinearity f t satisfies (|23)l . ([^b]) for all 

te [0,1]. 

Since any solution to (|3.ip is bounded above by 1 we may redefine ft(u) as a 
constant for u > 1. Thus we may assume that for some constant M > we have 

|/t(u)|<M Vuel, V*G[0,1]. 

We define a nonlinear map T e>t : X — > X as follows. Let (it,/?) € X. Then solve 

(3.2) Av = /3 2 v/ t (M + (9) inBi, u(l) + u'(l) = l on dB x . 

This problem has a unique solution which can be found for example by minimizing 

1 / /i iO «o « / 1 



(|V«| 2 +/3 2 / t (« + 0) v 2 ) +- / w 2 -/ v 

7 L -IdBi JdBx 
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(3.3) 



with v £ H 1 (Bi). Since ft(u + 9) e L°°(Bi) by standard elliptic regularity v € 
C 1,a ([0, 1]). Then find /3 > and u such that 

( -Au + e(3 2 u = /3 2 vf t (u + 9)- e0 2 9 + /3 2 - /3 2 in Si 

\ u(l)=0, u'(l) = -fl(l + j8Ve) 

admits a solution. This problem has a unique solution (u, (3) which furthermore is 
C 1,Q ([0, 1]) (this assertion is verified in the proof of Lemma \3 . 1 1 below) . We define 

T £!t (uJ) = (u,/3). 

Observe that (it, /3) is a fixed point of T e t if and only if (u + 9, v) is a solution of 
©Ul- 
Leninia 3.1. Write T(u,P,e,t)=T e>t (u,P). Then 

T:lx[0,co)x [0,1] -> X 

zs compact. Moreover 

(3.4) ||T(n,/?,e,t)|| x < C(l + sinh(Vi/?) + /3 2 ) V(u,/3) e X.Ve > 0,Vt e [0,1]. 
Proof. First we remark that the solution w to p.2[) satisfies 

< u < 1 in Si 

because is a subsolution and 1 is a supersolution. Given a > let Z a be the 
solution to 

-AZ a + aZ a = infli, Z„(l) = 1. 
Note that Z a is explicit: 

„ , . sinh(ar) 
Z a (r) = . ).[ V0<r<l 
r sinn(o) 

and hence 

(3.5) — ^— < Z a < 1. 

smn(a) 

Let (u,(3) = T(u, /3,e,t). Then choosing a = s/3 2 , multiplying p.3[) by Z a and 
integrating in £?i we find 

-u'(l) = (3 2 f {vf{u + 9) + l)-p 2 f (e0 + l)Z a . 



^From the boundary condition in (|3.3[) it follows that j3 is explicitly given by 

6(l + y-eP) + p 2 J Bi (e9 + l)Z a 
( ' ' J Bi (vf(u + 8) + l)Z a ' 

where a = sj3 2 . By (|3~5|) 



(3.7) / („/(fi + 0) + l)Z a > , ff ) 

Jb 1 smh(y/ep 

^From formulas (|3.6|) . (|3.7[) we see that 

f3 2 < C(l + sinh(Vi^)+/3 2 ). 

Using this inequality, standard elliptic estimates and the facts that v < 1 and /* 
is uniformly bounded we obtain (|3.4j) . From here we deduce that T is continuous 
and compact. Indeed, for the latter assertion, note that if B be a bounded set of 
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X x [0, oo) x [0, 1] then T{B) is a bounded set in C^QO, 1]) x K for any fx g (0, 1) 
and taking /i > a it follows that this set is precompact in C 1,Q ([0, 1]) xR. □ 

Lemma 3.2. The operator T^o has a unique fixed point («o,/3o) WJ Af. 
Proof. In this situation u + v = 1 and hence the system reduces to 

-Ait = I3 2 {1 - u)(u - 0)+ in Si, u|e Bl = 
with the additional requirement that 

Let w = u — 9. Then the equation for w becomes 

(3.8) -Aw = /3 2 (l-0-w)u; + in B u w\ 9Bl = 0. 

This equation is of logistic type and many properties are well known (see [131120] ): 

a) Any solution w to (|3.8[) satishes < w < 1 — 0, and either w e or w > 
in Bi. 

b) ()3.8[) has a nontrivial solution if and only if j3 > (3* where /3* is such that 
(/3*) 2 (1 — 0) = Ai and Ai denotes the first Dirichlet eigenvalue of —A in 
the unit ball. 

c) For p > f3* there is a unique non-trivial solution, which we write as wp. 

d) uifj is monotone increasing with respect to /3 and 

lim wr = 1 — 

/3->oo 

uniformly on compact sets of B\ and 

lim wn = 0. 

13^/3* 

^From the above properties it follows that 

r dw P a v dw P n 
lim — = oo and lim — — = U. 

p^foo OV /3— 5-0 OV 

Moreover 

9wb . . 
p i— > — - — is strictly increasing. 
ov 

It follows that there is a unique (3 > /3* such that 

dW l3 = Q 

dv 

We call this value flo and the let uo — wp . Then (uq, /?o) € AT is the fixed point of 

r 0l0 . □ 

Lemma 3.3. Let (uq,(3q) denote the unique fixed point of Xo ; o found in Lemma \3.2i 
Then (uo,/3o) * s nondegenerate. 

Proof. Let us write (uo,0o) & X the solution of / — Xb,o = of Lemma T3. 21 We 
have to verify that the linearization of / — Xo.o around this solution is an invertible 
operator. Note that the operator X^o involves only the nonlincarity /o(s) = (s — 0) + 
and that the fixed point (uo,f3o) from Lemma 13.21 satisfies uq > in [0,1) and 
u' (l) < 0. Therefore for ip in a neighborhood of uq in the topology of C 1,a ([0, 1]) 
we have fo(uQ + ip + d) = (uo + ip) + — ug + ip. Hence Tq.q is differentiable at (uq,/3o)- 
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We next compute the derivative of Ti^o at (uq, Po) in the direction of {ip, a) £ X, 
which we write as 

WS7) = DT 0fi {u ,(3 )(ip,a). 

Then 

(3.9) ^ = ■02 + 0^1, 7 = 72 + ct7i 

where (V'2,72) and (-01, 71) are computed as follows. 

To compute ip2 we have linearize (|3.3[) with respect to u and then set u = u = uq, 
P = P = (3q. Since for t — 0, we have fo{uo + 0) = u o and Jq(uq + 6) = 1, the 
function ^2 satisfies 

^ 3 1Q v f - A^2 = 2f3 "f2VQU + Po^ u o + Po v oP + 2/?o72 hi B x 

\^(1) = 0, v 2 (i) = o, 

where 72 is adjusted so that -02 satisfies both boundary conditions, is found 
solving 



(3.11) 



9u 9r 9{t 



and t>o is the solution of (|3.2[) with /3 = /3o, u = ito and £ = 0. As explained in the 
proof of Lemma 13.21 then 

v = 1 - - Uq. 

To compute ipi we linearize (|3.3[) with respect to /3 and then set u = u = uq, 
P — (3 — /3q. Therefore satisfies 

dv 

AtPi = 2p 0ll v u + Pq-~u + 2/3 7i - 2^ in B x 
(3.12) { dp 

Vi(l) = 0, ^(1) =0, 

isen so thi 
puted from the equation 

(3.13) 

We need to find the kernel of I — D- zTofi(uQ, fig), that is solutions (<p,a) G X 

of: 

(3-14) £L jTq.oK *) = {<P, °) 

Using the notation above we may write ip = tp2 + o"ipi , & = 72 + 0*71 ■ We claim that 
m ir\ /dv dv s 

3.15 lf=- — r 

V9m flfl. 



where 71 is chosen so that both boundary conditions are satisfied and |~ is com 

dv dv 
A— -~ = 2p v u + Po—~u in Bi 
dp dp 

dp drdp 



dp 
cal 

Az = in Bi 



To prove this, let z = f + || + 17 Then a calculation shows that 
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and 2(1) + z'(l) = 0. This implies z = and (|3.15[) follows. ^From the formulas 
(|3l ^ - (prT3l) we get 

-Aip = Pl(l-e- 2u ) p + 2crft(l -9- u )u 

But also (p satisfies 

The linear elliptic operator — A— ft 2 (1 — 6 — 2uq) with Dirichlct boundary condition 
arises as the linearization of (|3.8[) around the solution Uq, and it is well known that 
it has a positive first eigenvalue. Since ft(l — 8 — Uo)uq > 0, for positive a we 
would have p > and ip' (1) < and for negative a we would have <p < and 
(p'(l) > 0, which are both impossible since p>'(l) — 0. Hence a = and also p — 0. 
This shows that the linear equation D- pTo t o(uo, ft)(v?, a) = (tp, a) admits only the 
trivial solution. Since ^Tq^uq, ft) is compact we see that I — pTopfao, /3q) 
is an isomorphism. 

□ 

Proof of Theorem II. li By Lemma \2.2\ there exists <5 > and /3i > such that 
for all (3 > P\, fi^fs < S and all t £ [0, 1] there is no solution in S to the system 
(|2.1[) . This means that T e t has no fixed point (u,/3) if fi > $\ and fiyfe < S. Let 

i?! = (7(1 + sinh(<5) + P'f) 

and define 

fit = {(u,(3) G X : ||«|| OI .a ([0il]) < Ri, <y9 < 

Let £i = {8/ Pi) 2 . Then f2i is a bounded open set of X and for < e < Si and 
< t < 1 the operator T e t has no fixed point in dfli. Indeed, suppose (it, 0) £ dfli 
is a fixed point of T £ t . It is not possible that (3 = Pi by Lemma [2.21 The case 
ft = is also impossible. This means that < /? < fi\ and ||u|| c i,c,([ ^ = R\. But 
by inequality p.4p we would have 

Ri = ||u|| o i,- ([0) i]) < C(l + sinh(Vi/3) + /3 2 ) < C(l + sinh(<5) + ft 2 ) = flj 

which is impossible. Thus the Leray-Schauder degree deg(I — T £it ,fii,0) is well 
defined for < e < ei, t £ [0,1]. By Lemmas 13.21 and 13.31 we have 

de5(7-To, 0i ni,0) = ±l 

and by homotopy invariance 

deg(I-T Sit ,toi,0) = ±l V0 < e < e u Vie [0,1]. 

This shows T e t has at least one fixed point in fli for < e < £1 and t £ [0, 1]. In 
particular the system (|2.1[) has a solution with bounded /3 for any < e < e\. 

We know there exists e* > such that for e > e* the system (|2.1[) has no solution, 
and hence T e> i has no fixed points in X for such e. Let < £0 < £1 an d define 

«2 = {(u,P) £ X : [|ix|| cl , a([0)1]) < R 2 , < P < P2} 

where 

R 2 = (7(1 + sinh(V2F/3 2 ) + /?|) and ft = 6»/Vi . 
Then T e j has no fixed points on dVl2 for £o < £ < 2e* and i S [0, 1]. Then using a 
homotopy along e £ [eg, 2e*\ we find 

deg(I -T £A ,Q 2 ,0) = Ve <e<2e*. 
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This implies that for eo < £ < £1 the system has at least another solution and that 
this solution has /3 in the range 5/y/e < (3 < 9/ \fe. Since £o is arbitrary we obtain 
the same conclusion for < e < e\. □ 

4. Existence when / is discontinuous 



In this section we prove Theorem 11.21 We assume that / satisfies (|1.2j) . (|1.4j) . 
(jO]) and that 

j] = lim f(u) > 0. 

For n > 1 define 

f n (t) = mm(f{t),n(t-6)+). 
Then /„ is continuous and satisfies (|1.2[) and () 1 .3[) with a fixed constant. Moreover 
there is ho satisfying (|2.2[) and (|2.3[) and such that /„ > ft-o for all n. Such ho can 
be taken for instance as 

h (t) = min((f - 9) + ,a) 

for some a > small enough. 

By Theorem 11.11 there is e% > such that the system (|2.1I) admits 2 solutions 
(it*,^), (m^,^) for all e € (0,£i). Moreover there is a fixed number /?i such that 
Pn < A and /3 2 € (S/y/e, 1/(9^/5)) for all n > 1 and ee (0,£i). From now on we 
fix e € (0,£i) and study the limit as n — > oo of any of the 2 solutions which we 
call just (u n ,v n ) with parameter /3 n . Since /3„ is bounded and by standard elliptic 
estimates we may assume that /3„ — >■ ft, u n — > u and v n — > v in C 1,Q ([0, 1]). 

Since f n (u n ) is uniformly bounded we may extract a further subsequence such 
that f n {un) — > <P weakly-* in i°°(0,l). The function tp then satisfies < <p < 
max[ 01 ] / in [0, 1]. Since u n satisfies a linear equation for r > 1 and vanishes at 
infinity it has an explicit formula 

e -(r-i)£„Vi 
u n (r) — 9 for all r > 1. 

Setting ip(r) = for r > 1 we see that u, v satisfy 

Au - £/3 2 u = -(3 2 vip in M 3 
Au = /3 2 v^ in M 3 
lim = 0, lim v(x) = 1 

^oo \ x \— 

The set D = {r£ [0, oo) : u(r) ^ 9} is open (relative to [0, oo)) and for any r E D 
we have f n (un(f)) — ^ f( u ( r )) asm oo. It follows that ip = f(u) in D. Outside D 
we have (p < r\ a.e. Indeed, let ip £ C^°((0, 1)), ip > 0. Then 

fn(Un)tp= / fn(u„)l(>+ / f n (u n )lp. 

JDn[o,i] JD c n[04] 
Since f n (u n ) — > /(m) in I? by dominated convergence we have 

fn(u n )lp -> / 

'Dn[0,i] JDn[o,i] 

On the other hand, if r G Z? c limsup„_ >+00 f n {u n (r)) < r) so that {fn{u n )^f]) + — > 
on L» c . Then 

/ (/«(««)- I?)^ = / (fn(u n )-r)) + i;- {fn{Un)-Tl)~Tp < 

JD c n[o,i] Jfl c n[o,i] Jr» c n[o,i] 
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where o(l) denotes a sequence converging to as n — > +00. It follows that 

fn(u n )4> < rjTp + o(l) 

'L> c n[o,i] JD c n[o,i] 
and hence 

tpip = lim / f n {u n ) < I f(v,)ip+ / 

la n ^+°°Jo JDn[o,i] JD=n(o,i] 

This shows that ip < r\ a.e. in D c . 

Our main objective is to show that the complement of D is finite. If u'(r) 7^ 
whenever u(r) = 9 then D c is discrete and since it is contained in [0, 1] it is finite. 

Let us analyze the case where for some re [0, 1] we have u(r) = 9 and u'(r) = 0. 
Let 

r = sup {r > : u(r) = 9, u 1 (r) = 0}. 

Then ro < 1 and u(ro) = 9, u'(ro) = 0. We assert that there is a small interval 
(^Oj r o + c), c > such that u > 9 in (ro,ro + a). To prove this we start ruling 
out the possibility that u{r n ) = 9 for some infinite sequence r n \ ro. We actually 
may assume that if n is even then u > 9 on (r n+ i,r„) and if n is odd u < 9 on 
(r„ + i,r„). Let us see that the following holds 

(4.1) u'(r n+1 ) = -u'(r n ) + 0(£ 2 n ) 
where 

and 0(£ n ) denotes a sequence bounded by Cl 2 n with C independent ofnasn^oo. 
Suppose first that u < 9 on (r n+ i, r n ) and define u n by 

u(r) = fl+^u n ((r-r n+1 )/^ n ) re [r„ + i,r„]. 

Then 

op 

< + -T— ~ < = £ P ( £ nUn + 0) for r e [0, 1]. 

u n (0) = u n {l) = 0. 

This equation implies that u n and u' n are uniformly bounded on [0, 1]. Therefore 
< = ef3 2 9 + 0(£ n ) for re [0,1] 

hence integrating 

«n(r) = ^£/? 2 (l + 0(^ n ))r(r - 1) for r e [0, 1] 

and 

(4.2) <(r) = ie/3 2 (l + 0(£„))(2r-l) for re [0,1]. 



This proves (|4.1[) in this case. 

Now suppose that u > # on (r„ + i, r n ). Defining u n as before we now obtain the 
following equation 



u'„ = ef3 2 {£ 2 n u n + 9)- P 2 i n f{9 + £ 2 n u n ) for r e [0, 1] 



n r^„ + r„ +1 

«n(0) = G„(l) = 

where 

v„{r) = v(r n+1 + t n r) r e [0, 1] 
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This equation implies that u n , u' n are uniformly bounded on [0, 1]. Thus 

< = e/3 2 6 - {3 2 v n (0)f(9 + i 2 n u n ) + 0{l n ) for r e [0, 1]. 
Multiplying this equation by u n and integrating in [0, 1] yields 

which shows that 

fin(l)' = -U n (0)' + 0{i n ) 

and proves (|4.1[) in this case. 
Using (|4.1[) inductively we find 

n+k—l 

u'(r n ) = (-l) k u'(r n+k )+ Vn>l, fc>l. 

j=n 

Letting k — > +oo, and using that u'(r ) = yields 

oo 

u'(r„) = ^0(£ 2 ) Vn>l. 
Choose a subsequence rij — > oo as i — > oo such that £ n; > £j for all j > rij. Then 

oo 

|w'(r rai )| < ^ £ J- 
But flO) implies that |u'(r„)| = £ n ej3 2 {l + 0(£ n ))/2. It follows that 

OO 

£ ni ep 2 {l + 0{t ni ))/2<C£ ni Y,£i 

which is impossible as i — > oo. 

This establishes the assertion that u cannot oscillate infinitely many times around 
9 to the right of tq. It follows that for some a > either u > 6 in (?"o, ?'o + c) or 
u < 9 in (ro, ro + er). The latter can in fact not occur by the Hopf lemma. 

Since u satisfies the ODE u" + |it' = e(3 2 u — (3 2 vf(u) in (ro, ro + a) and in this 
interval u > 9 we see that u is C 2 in [ro,ro + er). Since ro is a minimum of u 
restricted to [ro, ro + a) it follows that lim ? ,_ >j ,+ u" > 0. Which yields the following 
inequality 

(4.3) el3 2 9 - P 2 v{r )ji > 0. 

Now, if v is constant in [0,ro] this means that ip — a.c. in [0,ro] and u < 9 in 
[0,ro]. Actually u < 9 in [0, ro) by the strong maximum principle and this finishes 
the proof in this case. If v is not constant, then v(r) < v(ro) for all r € [0,ro). Let 
us verify that D c is discrete. Suppose that for some r\ € [0,ro) we have u{r\) = 9. 
If u'(ri) then r\ is isolated. If v!(r\) = then integrating the equation for u 
we have 

r 2 u'{r)=( t 2 (eP 2 u- /3 2 v<p)dt. 
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For r near n we have 

if u(r) < 9 then ef3 2 u(r) - f3 2 v(r)ip(r) > 

if u(r) = 9 then e/3 2 u(r) - (3 2 v(r)ip(r) > e(3 2 9 - /3 2 v(r )r] > 

if u(r) > 9 then e(3 2 u(r) - f3 2 v(r)ip(r) = e(3 2 u(r) - (3 2 v(r)f(u(r)). 

In the second line above we have used that v(r) < v(ro), if < r\ and (|4.3|) . In the 
third line above we may say that e/3 2 u(r) — (3 2 v(r)f(u(r)) > if r is sufficiently 
close to r\ by continuity, v(r) < v(ro) and (|4.3[) (we may regard / as continuous 
here since we are working with values above 9). This shows that u is strictly convex 
in a neighborhood of r\ and hence there are no other points in D c close to r\ . This 
shows that D c is discrete, hence finite, and finishes the proof of the theorem. 



5. The Heaviside ignition 

In this section we perform explicit calculations for the ignition nonlinearity 
f( u ) = X[u>e]i where < 9 < 1. We first reduce the differential equations to 
a finite number of equations in some parameters. In a second part we obtain rigor- 
ously the bifurcation diagram of Figure [T] for small e. 

We rewrite (12.11) in the form 



(5.1) 



and 



(5.2) 



Au- f3 2 eu = in(l,+oo) 
Av = in(l,+oo) 

u(i) = e, «(i) = 7 

m(+oo) = u(+oo) = 1 

Am - P 2 eu = -P 2 vf(u) in (0, 1) 
Av = f3 2 vf(u) in (0, 1) 

u(i) = e w(i)=7 

u'(0)=0 u'(0)=0 

u'(l-) = u'(l+) v'(l~) = v'(l+) 



where 7 is a parameter to be adjusted. 

Observe that problem (|5.1[) can explicitly be solved. Indeed, the following func- 
tions solve (15.111 



Hence the last condition in 

«'(!+) : 



r 

becomes 



v = l- 



1-7 



-B{l + Py/e) w'(l+) = l-7. 



5.1. Solving (|5.2j) with the assumption u > 9 in B(Q, 1). Assuming that u > 9, 
in -B(0, 1) the subsystem (|5.2[) is reduced to the following linear problem 

P 2 eu = -p 2 v inS(0,l) 



(5.3) 
(5.4) 
(5.5) 
(5.6) 
(5.7) 



Au 

Av = fi 2 v in B(0, 1) 
u{\) = 6 w(l)=7 

= 1- 7. 
u'(l-) = -e(l + /9Ve) 
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For e 7^ 1, observe that the two following functions 
1 sh(pr) 



(5.8) 



(5.9) 



pch(p) 



re (0,1) 



(a i th (0) \ 



re (0,1) 



1 - e sh{Py/e) r 
solve equations (|5.3l) - (|5.6j) . ^From the formula for v we can see that 

th(f3) 



P 



7 



Now to obtain a solution of subsystem (|5 . 2[) . it then remains to adjust j3 and e to 
get (|5.7p . This equation is 



-9(l+PV^)=u'(l-) = ^- 

1 — c 



/3(l-e) 



which after simplihcation yields 

(5.10) 6(l-e)(l + th(Py/e)) 
Let us define 

(5.11) g(e, p) := 0(1 - e)(l + i/i(/3^)) + 



th{p) th(fiy/e) 



P PV^ ■ 



P 



PV~e 



It can be shown that there is e* > such that for e € (0,e*), g(s, ■) has 2 zeros, 
f3~ < fif . The function <j(0,-) has a unique zero /3q = lim e _yoAr- A calculation 
(see Lemma IBTTj) shows that for e > small and P = p~ the solution it obtained by 
the formula (|5.9j) satisfies it > 6* in [0, 1). Along the lower branch (3~ , < e < e* 
we observe numerically that it > 9 in [0, 1). After reaching e* we continue the curve 
along the upper branch /3j now decreasing e from e* to until we reach a critical 
value e = £o where the corresponding solution given by (|5.9|) satisfies u > 9 in (0, 1) 
but it(0) = 9. 

5.2. Solving the subsystem (|5.2[) with the assumption it > 6* in 5(0, 1) \ 
5(0,7?) and u < in 5(0,7/) for some 77 > 0. 

Under this assumption the problem (|5.2[) decouples again into two subsystems 

/3 2 eu = -/3 2 v 
(i 2 v 



(5.12) 



and 



Ait 
Av - 
u(l) = 9 = u(rf) 
u'(l-) =u'(l+) 
u'(r)-) = u'( V + ) 



in 5(0, 1)\ B(0, rj) 
in 5(0, 1)\ B(0, rj) 

7j(l) = 7 1>(?/) = (5 

w'(l-) =v'(l+) 
v'(r ] -)=v'(r ) +) 







in S(0,?7) 
in 5(0,7?) 



!Au - /3 2 eit 
Av = 
11(77) = 77(77) = <5 

where S, 7, 77 are parameters to be found. Using the conditions at r — r\ and r = 1 
we will be able to reduce the parameters to only /3 and 77 which will be implicitly 
defined as functions of £ by 2 equations. 
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Observe that since v is harmonic in B(0, 77) with the Dirichlet condition v(r]) = 5, 
we deduce that v = 6 and therefore v'(r)—) = 0. 

Let us now solve (|5. 12p explicitly. First, observe that the function v = A sh ^ r ^ 4- 



D 



ch(Pr) 



satisfies 



Av-fi 2 v = in B(0,l)\B(0,rj). 
Therefore, choosing A and B such that 

1 = Ash(j3) +Bch((3) 

1 - 7 = A[/3ch{/3) - sh(P)} + B[f3sh(/3) - ch(/3)], 

it follows that v satisfies the right boundary conditions on dB(0, 1). Solving for A 
and B wc find 

A= and B = . 

To obtain a full solution it remains to adjust 7 in such way that v'(j]) = 0. But 
, A[ Pch(Pr) shjpr) or gfW gMM i 

V = A t " \+ B [ J 

and setting v' (rf) = in the above equation yields 



/3ch(/3r)) sh(Prf) 



V 



>r 



B 



f3sh((3ri) ch(Prf) 



(r)Pch{0r)) - sh{f3ri)^ + B^rfsh^) - ch(0rj) 



= 
= 



Substituting now the values of A and B in the above equation, it follows 
(ch(l3) - jPsh(p)) C - (sh(p) - lPch{P)) D = 0, 



where 



Therefore, we have 
(5.14) 



C := ( ril3ch(Pri) - sh(/3r]) 
( r)l3sh(Pri) - ch{Prj) 

C - Dth(/3) 



7 



p(cth{p) 



D 



which after simplification reads 
7 = l(P,v) 



i](3 - th(?]P) - th(/3) 


r)pth(Prj) - 1 




P 


(pp-th{riP))th(P) 


+ 1 - r]Pth{r]P) 



Thus we may compute now S = v(rj): 

' ch{P) - jPsh(p) \ sh{/3r]) ( sh(f3) - jpch(P) \ ch{prj) 



5 = 



P 



V 



V 



^From (|5 . 1 3[) we know that in B(0,rj), we have the following 

6ri shirpJe) 
u = . 
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Therefore, 

" /3^ich(r)/3^I) 1 



(5.15) u'{rj) 



sh(r]/3y/£) 7? 



rjsh(?]/3 v 'e) 



Let us now find the solution u in B(0, 1) \ i?(0, 77). As for the construction of v, 
let us observe that for e 7^ 1, 

(5.16) U: =^ +£ ^M) + F^M) 

e — 1 r r 

solves the following equation 

Au - f3 2 eu = -f3 2 v 

By taking 

(5.17) E(fi,r,,i) := _ ,( c , (/3Vi) + s , ( ^)) 

(5.18) F(/3, 77, 7) := _fMW|_7|^(M) + e{ch { ^) + .AflVi)) 

we easily verify that u(l) = 6* and u'(l _ ) = — 0(1 + \/e0)- It remains now to impose 
the boundary condition on <9-B(0, 77). From the formula for u (|5 . 16[) we have 

uvn) = + h h r 

£ — 1 77 77 

(5.19) «'(„) = EfjyTe^^ + Ffjy^'W^ E sh ^f ll] - F cM/? f ^ 

77 7/ 7/ 2 77 2 

since 7/(77) — 0- Thus the equation 77(7/) = becomes 

(5.20) -*- + E(fi, 77, 703, + F(/3, 77, 7 (A 77))*^!) = , 
e — 1 ?/ 77 

and (|5.15[) combined with f|5 . 19[) becomes 

E(f3, 7/, j(f3, 77))77/3Vi c /7(/3yi77) + F(/3, 77, 7 (/3, 77))t7/3^ s /i(/3V^7 7 ) 
-£(/3, t?, 7 (/3, ri))ah{fiy/hi) - F(f3, 77, 7 (/3, rf))ch{^rf) 

1" 



(5.21) 



We solve numerically equations (|5.20[) and (|5.21j) with /3 and 77 as unknowns 
that depend on e, for < e < £o, where £0 is the critical value of e described in 
Subsection 15. II The result from this numerical computation is shown in Figure [TJ 

5.3. Results for small e. 

Lemma 5.1. Fore > suitably small /3~ is a smooth function of e and the solution 
given by (|5.9[) corresponding to (3 — /3~ satisfies u > 6 in [0, 1). 

Proof. Let fti the unique solution of the equation 

thjp) _ (1 - g) 
[3 2 ■ 
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Then 



g(e,/3 1 ,6)=6(l-e)(l + th(0iV^)) + 



1-6 _ tMfixy/e) 
1 + 6 -26s th(P ly /e) 



= 6{l-e)th{p^)) 



1 - 



-6e + [l- 



2 ftVe 
For p ly /e < 1, using that fiiy/e - (/?1 ^ )3 < th(fi ly fe) < fiiy/e, we end up with 

9(e, 0i,O) < 6(1 e)PiVS- ^ - 6s + 



< 



(/3iv^) 2 



1 - 



Computing now the positive roots of + 6*X — yields g{e,P\,6) < for 



< 

Therefore, we achieve g{e,px,6) < 0, for e in (0, min{ 
Let us denote 



30 



/i i 2(1-6 



2/3i 



El 



30 a/1 



2(1- 



30 2 



1^ 



V 



2/3i 



/ 



Hence, for any (e, 6) in (0,min{ei, 4z}) x (0, 1) fixed there exists two possible so- 
lution to g(e,j3,6) = 0. Moreover, /3~ < jSi < /?+. 

To obtain a solution, to problem (|5.2[) , we still need to show that the constructed 
solutions effectively satisfies the conditions u > 9 in £?(0,1). Let us observe that 
u — 9 satisfies: 

A(u -9)- sf3 2 (u -9) = -f3 2 {v - eff). 
If v(0) — e9 > then using the maximum principle, it follows that (u — 9) > 0. 
Since v(r) = ( jj^p^j j sfe (ft r ) , it follows that u(0) = ^gy- Hence, we end up with 
the condition 

9 

> £0. 



ch{P) 



We can conclude with a final estimate on e, namely since (3 < /3i, and is a 

decreasing function, we have the following uniform estimates on e 

r< (l-6)Pi 
~ 6sh{p{) ' 

We therefore have construct a solution to the problem (1), when (e, 6) belongs to 
(0,e (6>)) x (0,1) where 



EXISTENCE OF FLAME BALLS 



25 



□ 

Lemma 5.2. Fore > small the system (|5.20j) . (|5.21j) has a unique solution f3j ,r) 
with i] > 0, fif is a smooth function of e and u defined by (|5.16j) satisfies u > 8 in 

Proof. The 2 equations f|5 . 20f) and (|5.21j) can be written in the form 

ijd 



- £ 
a i] 



tanh(a?7) 



Esmh(an) + F cosh(an) = r/6 + - 

E (arjcosh(an) — smh(ar])) + F (arj smh(ar)) — cosh(a?7)) = r\0 
where 

a = Py/e, 

Instead of e, f3, r\ consider the variables t, a, x defined by the following relations 

x 

t = y/e, a = fiy/e, r) = 1 - -. 

Consider also the function G(t,a,x) = (Gi(t,a,x),G2(t,a,x)) where 

r)S 



Gi(i,a,x) = 



Esinh(ar)) + F cosh(an) — r)6 



1 

G2{t, a, x) = E (an cosh(a?7) — sinh(a?7)) + F (an sinh(a?7) — cosh(a^)) 
\tanh(a77j 

where 5, 7, E, F are as before. We see that we have a solution to the system if and 
only if G(t, x, a) = 0. 
Using the fact that 

lim P'y = t~tt an d lim f35 



t-»-0 tanh(x) t-to sinh(x) 

it is not difficult to verify that G maybe extended in a C 1 manner for t = 0, 
a, x € (0, 00) with the values 

G\ (0, a, x) = —x + Oax + 



tanh(x) sinh(x) 
G 2 (0,a,x) = l-6(l + a)-9(^-l 

Note that the equation G 2 (0, a, x) = is equivalent to 

tanh(a) J- 

1-61(1 + a) ~~ 6 

which is seen to have a unique solution ao G (0, 1/9 — 1) since the left hand side 
defines a strictly increasing function of a which goes to zero as a — >• and goes to 
+00 as a /■ 1/6 - 1. 

Now the equation G±(0, ao, x) = can be written as 

(5.22) 1- 7 * + \, s =g(l + o )6(0,l). 

xtann(x) 1 sinn(i) 
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The above left hand side is strictly increasing with limits 1/2 as x — > and 1 as 
x — > oo. Using the inequality 

1 < a Va > 



tanh(a) 

we see that 1 — — 9a < 6a Q which yields ao > Hence 

1 6» 1 
0(l + a o )>- 2 + - 2 >- 2 . 

Thus (|5.22j) has a unique solution xq S (0,oo). □ 



6. Estimate on e* 

As mentioned before the largest value of e such that p. II) has a nontrivial solution 
satisfies 

e* <e° 

where 

e° = inf{e > : ff e (s) < Vs € [0, 1]}, 

and 

ff e (s) := (1 - s)/(s) - es. 

In the particular case of /(w) = X[u>e], < 9 < 1 we obtain the value e° = 4. This 
estimate can be sharpened. 

For this, given e > let us introduce 

G e (b)= [ g E (s)ds. 
Jo 

For < e < So define also a(e) as the smallest zero of g e in the interval (0, 1). Wc 
b(e) be the largest b in [0, 1] such that g e is positive on (a(e), b(e)). 

Proposition 6.1. Assume that f satisfies (|1.2[) and is continuous. Then (| 1 . 1 1) /ias 
no nontrivial solution if £ > e 1 where 

(6.1) e 1 = sup{e > | G e (b(e)) > 0}. 

Before proving the above proposition, let us recall a result on the scalar problem 
on a bounded domain il C K w : 




—Aw; = h(w) in fl 

w = on dil. 



There is a vast literature concerning the existence of positive solutions of the above 
equation. Wc will just mention a necessary and sufficient condition for the existence 
of positive solutions. 

Theorem 6.2. ( Clement and Sweers |14j j Let fi C M. N be a bounded smooth domain 
and h € C (R). Assume there exists &i < 62 with 02 > such that h{b\) = h{b%) = 
and h > on (61,62)- Then there exists positive solution w to (|6.2|) with maxiu £ 
(61, 62) */ o^rf onZy i/ J^ 2 h(s)ds > /or 6 G [0, 62)- 
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Proof of Proposition KTTl Fixe > e . From our assumption, we have G (b(e), e) < 

0. Let g a smooth function such that g £ > g e and J g £ (s) ds < 0. This is always 
possible since G(b(e),e) < 0. Observe now that from the necessary and sufficient 
condition given in theorem 16. 2\ the following scalar problem 

I Aw = —g e (w) in f2 
I w = on dfl 



(6.3) 



does not have any positive solution for any bounded domain f2. 

Now, we argue by contradiction, assume there exists a non trivial couple (u, v) 
solution of the system An easy computation shows that u is a sub-solution 

of the problem «E3) with SI = W l . Take 5 > small, and define 



9e,S =9e{s + S). 

So for S small enough, (i.e 5 < So), we have: 

rb{e)—5 rb(e) r& 

G s (b(e) -S)= j g £ ,s(s)ds = / g e {s)ds = G{b{e),e) - / g e (s)ds < 
Jo Js Jo 



Therefore g e ^s does not satisfy the necessary existence condition of Theorem 16.21 
Now, take Q — B(0,R). Since u goes uniformly to when \x\ goes to infinity, we 
have sup af2 it — > as R — > +oo. Thus, with 5 = sup^Q u < So, us = u — 5 is a 
sub-solution of the following scalar problem: 



(6.4) 



Aw = -g s (w) mB(0,R) 
w = Q on dB(0, R) 



Observe that the constant b(e) — 5 is a super-solution of (|6.4[) and us < b(e) — 5, 
then we can apply the monotone iterative scheme to obtain at least one positive 
solution which contradicts theorem 16.21 □ 



In the case /(it) = X[u>6] where < 9 < 1 with a similar argument we can show 
that for e > - — t>>— the system (jl.lj) has no nontrivial solution. Indeed, given any 
a > choose a smooth function f > f such that f(u) = 1 for u > 9 and f(u) = 
for u < 9 — a. With the same argument as before, (jl.lj) has no nontrivial solution 
for e > e 1 where e 1 is given by (|6.1[) with / replaced by /. A computation then 
shows that as a — » 0, e 1 — > ■ 
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